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SUMMARY 


A multivariable synthesis procedure based on eigenvalue/eigenvector assign- 
ment is reviewed and is employed to develop a systematic design procedure to 
meet the lateral handling qualities design objectives of a fighter aircraft 
over a wide range of flight conditions. The study reveals that the closed-loop 
modal characterization developed provides significant insight into the design 
process and plays a pivotal role in the synthesis of robust feedback systems. 

The simplicity of the synthesis algorithm yields an efficient computer-aided 
interactive design tool for flight control system synthesis. 


INTRODUCTION 

It has been recognized that for flight control system design, state space- 
oriented procedures provide significant design advantages over classical con- 
trol design using root locus. Bode, and Nyquist techniques. If the aircraft 
model is represented in the state space framework, two well-known synthesis pro- 
cedures are the linear quadratic regulator design and pole-placing technique. 

The quadratic synthesis procedure is particularly well suited to those systems 
which have desired response time histories as performance specification. How- 
ever, aircraft handling qualities design objectives originate in the frequency 
domain as desirable locations of closed-loop poles (eigenvalues in state vari-^ 
able representation) with satisfactory damping characteristics. 

A major difficulty encountered in the quadratic regulator formulation is 
selection of appropriate optimal performance indices to meet desired handling 
qualities. This difficulty has been partially resolved by characterizing the 
quadratic performance indices in terms of asymptotic eigenvalues and eigen- 
vectors. (See ref. 1.) This study (ref. 1) confirms that the modal control 
formulation is an appropriate synthesis framework for aircraft handling quali- 
ties design. 

The synthesis of multivariable systems based on pole-placement criteria 
has developed into several well-documented procedures (refs. 2 to 4) but none 
have proved to be effective design procedures. The common difficulty is not 
the fault of the algorithms but is really due to the inherent inability of pole 
(eigenvalue) specifications to characterize the actual variable responses. The 
state feedback law assigning a specified set of eigenvalues is not unique for 
multi-input systems. Different laws yield identical eigenvalues but radically 
different eigenvectors. Further, an eigenvector determines the influence of 
the corresponding eigenvalue on the state variable response. Thus, control of 
the modal matrix (matrix of eigenvectors) structure becomes as essential as 
satisfying eigenvalue specifications in order to effectively shape the dynamic 
response of the feedback system. 

Based on this key concept, it is possible to formulate the multivariable 
synthesis problem as an eigenvalue/eigenvector assignment problem. (See 



ref. 5.) This new formulation reveals that for an n-state, m-input system, up 
to m entries in each closed-loop eigenvector can be arbitrarily chosen, in 
addition to assigning all n eigenvalues. Thus, the eigenvector selection 
freedom can be successfully used for response shaping. 

The purpose of this report is to establish systematic synthesis procedures 
for the selection of eigenvalues and eigenvectors to meet a wide variety of per- 
formance requirements. In particular, the eigenvalue/eigenvector modification 
process is illustrated. Control laws for a typical fighter aircraft are synthe- 
sized to improve Dutch roll damping, turn coordination, and sideslip gust 
response over a wide range of flight conditions. 

SYMBOLS 

Values are given in SI and U.S. Units where appropriate. Any consistent 
set of units may be used elsewhere. 

A system matrix 

a, ai scalars defined in equations (Al) and (A5) 

3 y lateral acceleration, m/sec^ (ft/sec^) 

B input matrix 

b, bj^ scalars defined in equations (Al) and (A5) 

C, C measurement matrices defined in equations (29) and (Bl) 

Ci,Cj^ matrices defined in equations (10) and (17) 

D, D matrices defined in equations (29) and (Bl) 

dij scalar defined after equation (23) 

er vector defined in notation (1) of "Spectral Synthesis Algorithm" 

F matrix defined in equation (8) 

fj. vector defined in step 4(a) of algorithm 

G matrix defined in equation (8) 

g acceleration due to gravity, 9.8 m/sec^ (32.2 ft/sec^) 

g^(k) vector defined in step 4(a) of algorithm 

H matrix defined in equation (33) 

hr vector defined in step 4(a) of algorithm 
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^index (index) th order identity matrix 

i,j indices 

/s ~ 

K,K,K feedback matrices 



feedback from jth variable to aileron channel, j = p,4>,B,ay,4) 
feedback from jth variable to rudder channel, j = p,4>,8,ay,(j) 
aileron to rudder interconnect gain 

elementary permutation matrix defined in observation 4 following 
algorithm description 


L^,Lr,Lg,Lg^,L 5 ^ aircraft derivatives defined after equation (29) 


Mp,Mj.,M 3 ,Hj),M 5 ^,M 5 ^ aircraft derivatives defined after equation (29) 
m(I^) matrix defined in equation (16) 


m number of inputs 

(k) (k) 

m^ ,mb vectors defined after equation (16) 

N matrix of eigenvectors defined prior to equation (15) 

Np,Nr,Ng,N 5 ^,N 5 ^ aircraft derivatives defined after equation (29) 
n number of states 


P matrix defined after equation (31) 

p roll rate, deg/sec 

Q(k-l) matrix defined in notation (6) of algorithm 

q vector defined in equation (Al) 

R(l^~i) matrix defined immediately before equation (15) 
r e (A(k)} r is an element of the set 

S matrix defined in equation (8) 

s vector defined in equation (Al) 

ti vector defined in equation (21) 

t time index 


u input vector defined in equation (1) 


3 



input vector defined in equation (2) 
modal matrix (matrix of eigenvectors) 


“P 
V 

Vq total equilibrium velocity, m/sec (ft/sec) 

Vi,vj,vjj eigenvectors 

vectors used in equations (8) and (A 5 ) ; a partition of v^jV^+i 
X state vector defined by equation ( 1 ) 

X G Rf' n X 1 real vector x 

Yp,Yr ,Y3,Y5^,Y5^ aircraft derivatives defined in equation ( 29 ) 
y,y output vectors defined after equations ( 29 ) and (Bl) , respectively 

^i' ^i+l vectors used in equations (8) and (A 5 ) ; a partition of v^jV^+i 
aQ^,a2,a3,ai scalars defined in numerical example 1 
3 sideslip, deg 

A(l),A(k) set of indices defined in notation ( 5 ) of "Spectral Synthesis 

Algorithm" 

6 scalar 

62^,62,53 matrix elements used in equation ( 14 ) 

6^ pilot aileron input 

<Sr pilot rudder input 

63 aileron angular deflection, deg 

Sf rudder angular deflection, deg 

6X|^,6z|^ perturbed quantities defined in step 2 and step 5 of algorithm 

e scalar 

£ 1 ,^ 2' ^3 matrix elements used in equation ( 14 ) 

C<3 Dutch roll mode damping factor 

9 sideslip parameter defined in quantitative aircraft performance 

specification 

A diagonal matrix of eigenvalues 
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^i'^k eigenvalues 

TTi mode condition number defined in equation (25) 

p 2x2 matrix defined in equation (A4) 

2x2 matrix defined in observation 1 following algorithm description 
Ojf scalar defined by notation (4) of algorithm 

roll subsidence time constant 
<t> bank angle, deg 

’ip yaw rate, deg/sec 

Dutch roll mode natural frequency 
Superscripts: 

T matrix transpose 

-1 matrix inverse 

transformed quantity 

Notations: 

II II matrix or vector norm 

II Euclidian norm 

I I absolute value 

» far greater than 

« approximately 

TT continued product used in observation 4 following algorithm 

Upper case letters of the alphabet denote matrices; matrix subscripts 
indicate partitioned quantities. A dot over a quantity denotes a derivative 
with respect to time. 

SPECTRAL CHARACTERIZATION 

The concept of rationally utilizing the freedom in control law selection 
to satisfy eigenvalue as well as eigenvector specifications forms a useful basis 
for the development of multivariable synthesis algorithms. Unfortunately, the 
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process of coupling the effect of the individual nodes (eigenvalues) to the 
state variable responses through the eigenvector entries is nonlinear. This 
fact becones apparent if the state variable response is derived by using the 
modal canonical form. However, it is possible to identify certain desired 
closed-loop system structures and consequently associated modal structures. 

For example, it is often the case that higher order systems may be considered 
as a coupling of lower order subsystems, each with its own specifications of 
acceptable performance. In such a case, the eigenvectors should be selected 
so that the eigenvalues appropriate to one set of response variables do not 
unduly influence the other responses. Similarly, it may be desirable to segre- 
gate short-time constant variables from long-time constant modes or to insure 
that systems with both real and complex pair poles have eigenvectors selected 
so that minimal oscillatory behavior will arise in those responses associated 
with real eigenvalues. 

The eigenvectors may, of course, be modified without disturbing pole loca- 
tions through the nonuniqueness of the modal control process. There is also 
insufficient freedom to completely select the individual eigenvector shapes. 
Thus, the general synthesis approach will be to structure the design process 
so that maximum capabilities are achieved for satisfying whatever eigenvalue/ 
eigenvector specifications exist. The multivariable synthesis problem will now 
be formulated as an eigenvalue/eigenvector selection problem as follows. 

Consider the controllable system 


X = Ax + Bu 


( 1 ) 


where x e Rn^ u e Rin, and A and B are properly dimensioned. Assume that 
B is full rank. The state variable feedback law takes the form 


u = Kx + Up (2) 

where K is the feedback matrix and Up is an external reference input. The 
problem is to select a K so that the closed-loop system matrix A + BK 
satisfies 


[a + BK]vi = XiVi 


(i = l,2,...,n) (3) 


where Xi is the ith eigenvalue and vi is the corresponding eigenvector. 

The formulation given by equation (3) permits a certain freedom in the 
choice of the eigenvector elements. In order to identify which of the eigen- 
vector elements can be arbitrarily chosen, partition equation (3) as 
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^ 


All 

Ai 2 


~Bi ■ 



~Zi 


Zi 

• • 

A2I 

• • 

A22 

+ 

• • • 

B 2 

[Ki : K2] 

> 

Wi 

= ^i 

• • • 

Wi 


'O- — -i j 


(i = 1,2,. ..,n) (4) 


where A^, B^, and are m x m matrices and other matrices are compatibly 

dimensioned. Also, assume that is nonsingular (if necessary by at most 

reordering the state variables), is an m x 1 vector, z^ and are 

partitions of the eigenvector v^, and v^*^ = [z^^ : . Completing the 

multiplication of the partitioned quantities yields 


[Aii + BiKi]zi + [Ai 2 + BlKaJwi = ^i=^i (5) 

[a 21 + B2Ki]zi + [a 22 + B2K2]wi = XiWi (6) 

From equation (5), it follows that 

KiZi + K2Wi = Bi"l[XiZi - AiiZi - Ai2Wi] (7) 


Substituting K^Zi + K 2 Wi from equation (7) into equation (6) yields the 
following fundamental eigenvector constraining relationship 


[^i^n-m - ^wi = [g + ^iS]zi 



»n) 


( 8 ) 


where In-m “ m) th order identity matrix and S, G, and F are 

matrices defined by 

S = B2Bx"1 

G = A 21 - SAii > (9) 


F = A 22 - SA 12 J 


Appendix A develops relations similar to equation (8) for complex pair eigen- 
value computations using real arithmetic. Now equation (8) constitutes a set 
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of n - m linear equations in the n unknown elements of each eigenvector. 
Thus, if is not an eigenvalue of F, at most m elements corresponding 

to the vector can be arbitrarily chosen. Then the remaining n - m eigen- 
vector elements corresponding to the vector can be computed from equa- 

tion (8) as 


Wi- = Ciz 




(i = 1,2, ...,n) (10) 


where = [Xiln-m “ ^i®l defined as the "modal coupling matrix" 

corresponding to X^. 

The preceding analysis shows that all n eigenvalues and up to n . m 
eigenvector elements can be arbitrarily assigned through the use of state vari- 
able feedback. The resulting modal matrix 


= [vi 


V2 


'n] 


however, must be nonsingular. 

It is also important to note that equation (8) shows the complete relation- 
ship between the eigenvalues and associated eigenvectors. This relationship 
not only illustrates the eigenvector forms that can be achieved by state feed- 
back, but it also points out immediately that certain preconceived objectives 
may be impossible to achieve. Later examples illustrate this point. 

Hence, the design approach is complete in the sense that when it fails to 
satisfy certain specifications, it does so by showing that no state feedback 
control law could satisfy them. Thus, the analysis provides a complete spectral 
characterization of all closed-loop realizations given the triple (A,B,Xi). 

It should also be noted that the free n . m elements of the modal 
matrix V (zi-vectors) correspond exactly to the n . m arbitrary elements 
of the feedback matrix K. This establishes the parametric equivalence between 
the nonunique feedback matrix K and the arbitrary modal entries. 

With a nonsingular modal matrix V and eigenvalue matrix A chosen 
subject to the constraints of equation (8), the closed-loop system matrix 
A = A + BK is uniquely determined by 


All • Ai2 


A = 


VAV'X 


A21 • A 22 


( 11 ) 
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The required feedback matrix K which yields this closed-loop matrix A 
can be easily computed by using the relations 


Ki = Bi~i[aii - Ai;[] 

K2 = Bi~i[ai2 - A12] 


(12) 


The preceding analysis relies on the generation of a nonsingular transfor- 
mation V (modal matrix). It is possible to prove the existence of a nonsing- 
ular modal matrix for any specified A (ref. 6). However, this proof imposes 
an unduly restricted structure on the modal matrix. This result is in contrast 
to a practical synthesis perspective where the designer may well be willing to 
accept a modest perturbation in a pole location if this will allow significant 
gains in flexibility of selection of eigenvector entries. This is particularly 
true if the designer will have control of the perturbation size and therefore 
will be able to control stability properties. This design philosophy has led 
to the development of the algorithm to be presented in the section entitled 
"Spectral Synthesis Algorithm." This algorithm consistently strives to satisfy 
both eigenvalue and eigenvector specifications while allowing the user to con- 
trol the perturbations in the design criteria. 

) 

The concepts of the modal control process will now be illustrated witji two 
numerical examples. The first example illustrates the eigenvalue/eigenvector 
selection procedure. The second example n»tivates the need for the development 
of a systematic algorithmic procedure to insure the generation of a nonsingular 
modal matrix. 


NUMERICAL EXAMPLES 


Example 1 


The following problem is presented to describe the synthesis procedure and 
indicates the trade-offs between design specifications and achievable results 
may be required. Suppose that the plant dynamics are 


*1 


1 

1 

-1 


*1 


1 

0 

_ 











ui 



0 

3 

-2 


X2 

+ 

0 

1 












U2 

X3 


1 

1 

-1 


*3 


0 

0 



. (13) 
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and the control law 


ui 

U2 


= K 


XI 

X2 

X3 


is to be selected so that 

( 1 ) Pole locations are at « -lOO, ^2 “ 3*^*3 X3 » - 1 . 

( 2 ) The variables X]^, X2, and X3 should exhibit, respectively, the 
short-, intermediate-, and long-time constant transients, namely, 0.01, 0.1, 
and 1 second, respectively. 

( 3 ) The responses should be decoupled in the sense that although X3 
may exhibit short and intermediate transients, e“^ will rapidly dominate 
other terms, X2 should not have e~^ terms, and xi should exhibit neither 
e“t nor e“X0t transients. 

Note that this is not a trivial problem since the open-loop system eigenvalues 
are at 0, 1, and 2 and the crucial long-time constant is to be associated 
with X3, the variable not directly influenced by the input. 

The time solution for the closed-loop system in terms of the modal 
matrix V can be expressed as 

x(t) = V exp(At) V“Xx( 0 ) 


or in terms of the individual components as 


Xl(t) 


vil V12 V 13 


aie^l^' 

X2(t) 

= 

V21 V22 V 23 


a2e^2t 

X3(t) 


_V 31 V32 V33 


036^ 


where the ot^ terms depend on V~X and the initial conditions x(0), and 
* - 100 , X2 « - 10 , and X3 » -1. Although the entries of the modal matrix 
influence the terms in a nonlinear fashion, dominance arguments show that 

the design specifications will be satisfied if V is of the form 
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^2 


V = 


Vll 

5l 


V22 

63 


^3 

V33 


(14) 


where the e and 6 values should be small compared with v^i (i = 1,2,3). 
This procedure will yield minimal cross coupling and preserve the diagonal 
dominance character of V“l. 


To begin the synthesis, denote each column of V by Vj, and 


Vj = 


^Ij 

^2j 

Wlj 


For the system represented by equation (13), the respective matrices of equa- 
tion (8) are 


S = 0 


F = -1 

G.[l 1 ] 


By requiring A «< Diag (-100,-10,-1) and starting with j = 3 (X3 «< -1 is 

the most crucial eigenvalue), the procedure is to select Z13, Z23, and X3 

and evaluate w^3 from 


W13 


Z13 + Z23 

X 2 + 1 


This is seen to be a "best case" problem (at least thus far) ; since X3 may 
be selected as exactly -1, both Z13 and Z23 selected to be zero, and 

wi3 =1. 

It should be noted that the eigenvalue selected coincides with the 
matrix F. This coincident eigenvalue selection will be discussed further in 
the spectral synthesis algorithm description. Thus, 
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0 


V3 



For j = 2, 


W12 


zi2 + Z22 
X 2 + 1 


where Z12 and W12 (corresponding to and 63 in eq. (14)) should be 

small compared with 222# and X2 should be about -10. For numerical sim- 
plicity, select X2 = -11, Z22 = ^12 ~ ^nd W12 is then found to 

be -0.1. Note that this selection of eigenvector entries is not "best case," 
since W]_2 snd 2^2 cannot be both 2ero if 222 is to be non2ero. This 
implies there must be some coupling of the system modes. By continuing the 
synthesis. 



It is seen by inspection that the two eigenvectors are independent. However, 
a higher dimension problem will require an algorithmic procedure to insure the 
linear independence of the eigenvectors. 

For j = 1, 


wii 


zil + Z21 


X + 1 


and for selecting X^ = -101, it is desired that zn dominate both w^ 
and Z21- Clearly, both conditions cannot be satisfied (again implying some 
coupling is unavoidable) . But by selecting = 1 and 221 = -1, w^ will 

be zero and the inverse of V will preserve the diagonally dominant structure. 
Again, 
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is clearly linearly independent of V 2 and V 3 by inspection. The modal 
matrix V has now been found to be 



The closed-loop system matrix is found to be 


A = VAV"1 



Thus, the required feedback gains are easily calculated from equations (12) as 
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The time response of the system is given by 


XI (t) 


~^ 10 ~ 

X2(t) 

= V exp [At]v^ 

*20 

_X3(t)_ 


_X3o_ 


Thus, 


XI (t) = xio e“101t 

X2(t) = -xio e"101t + (xj^Q + x2o)e"llt 

X3(t) = -O.Kxio + X2o)e~^^^ + (O.Ixiq + 0.1x20 + X3o)e“t 


It is seen that the performance specifications have been achieved and that 
the only coupling is from transients which rapidly become dominated by the 
desired time constant solutions. It is again emphasized that the example is 
"best case," since the ability to select four identically zero entries in the 
modal matrix insured the total unidirectional coupling. To now consider a 
"worst case" example, suppose the problem were the same except that 


A 



1 

3 

1 


-1 

-2 

19 


where only the element a33 in the A matrix has been changed. 

The first step of the solution is to again select X3 » -1, with 


wi3 


Z13 + Z23 
X 3 - 19 


and |wi3| » |zi3| and |wi3| » |z23|. These conditions certainly cannot 
even be approximately satisfied; thus, any solution will exhibit considerable 
mode coupling. 
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One particular compromise solution yielded a iiKsdal matrix 


V = 


-1 0 -10 

0-10 
0.05 0.033 0.1 


with A = Diag [j-1, -11, -8!^ , and the corresponding state variable responses 

XI (t) = (0.167x10 - 0.555x20 “ 16. 667x3o) 

+ 10(0.083x10 + 0.056X20 + 1- 667x30) 

X2(t) = X20 e“l^t 

X3(t) = 0. 05 (-0. 167x10 + 0.555x20 "*■ 16. 667x30) e“^ 

- 0.03x20 e“^^^ + 0.1(0.083x10 + 0.056x20 1. 667x30) 


The following observations may be made about this particular solution: 

(1) The solution for xi(t) is the only one which does not meet 
specifications . 

(2) If xi(t) were a more important variable than X2(t), selection of 
modal matrix entries could be made to have xi(t) free of cross-coupled 
transients (e“^,e~^®^), but this would introduce an e“^ transient in the 
X2(t) response. 

To summarize, the proposed technique yields desirable results in "best 
case" exan^les, and for "worst case" problems provides at the least a forewarn- 
ing of the type of mode couplings that will occur. Also, in less tractable 
problems the displaying of the system closed-loop modal structure shows that 
no other solution will yield the desirable response. These discussions also 
demonstrate the inevitable trade-off involved in any synthesis procedure and 
highlight the need to recognize fully the physical constraints of the plant, 
as exhibited by the modal coupling matrices (eq. (10)) to evolve accept- 

able design goals. 


Example 2 

A simple numerical example is now presented to highlight the final basis 
for the establishment of the design algorithm to be presented in the next 
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section. It will Illustrate that certain procedural idiosyncrasies, as simple 
as selecting the sequence of synthesizing eigenvectors, may induce problems 
which could be resolved by inspection for low order systems, but do require an 
algorithmic process for more complex systems. 

Suppose the plant dynamics are given by 


Xl~ 


~ 1 0 0~ 


XI 


~1 0 

*2 

= 

15 6 


X2 

+ 

0 1 

-*3_ 


-1 -1 0 


J'S. 


_0 0 _ 


“1 

U2 


with open- loop eigenvalues at 1, 2, and 3. Suppose further that design specif i 
cations are = -1, X2 = -2, X3 = -1, and an eigenvector structure 



-1 -1 

1 1 


wi W2 W3 


with the w elements computed to meet the desired pole specifications. 

To complete the design, equation (8) can be used immediately (where 
F = 0, S = 0, and G = [-I -3^ ) to yield 


1 
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It is immediately obvious that the resulting modal matrix V will be 
singular, since V2 and V3 are identical. For a third-order problem such 
as this, equation (8) may be examined directly for each eigenvalue, and a sound 
decision made as to which specif icat'ions should be relaxed and in what manner 
to give a system performance very close to that desired. Of course, for sub- 
stantially higher order systems or for an automated on-line controller, a pre- 
cisely defined algorithm is required. Such an algorithm is now presented, and 
this problem is reexamined as an illustration of its utility. 


SPECTRAL SYNTHESIS ALGORITHM 

As the discussions in the previous sections revealed, the central question 
in the new formulation of the pole-placement problem is the guaranteed genera- 
tion of the nonsingular matrix V satisfying equation (8) . An examination of 
the eigenvector constraints shows that there is an m-dimensional subspace asso- 
ciated with each eigenvalue; Thus, the problem reduces to selecting a nonsingu- 
lar set of n eigenvectors with each vector included from a subspace associated 
with a particular eigenvalue. However, this method is an inefficient way to 
synthesize since the designer loses direct control of the selection of arbitrary 
elements of the modal matrix (z-vectors). Also, any algorithm would become 
computationally intractable since it would involve pairing n vectors from a 
set of n . m vectors until a nonsingular set resulted. 

Alternatively, if the z-vectors are allowed to be chosen arbitrarily, then 
it is important to keep track of the linear independence of the eigenvectors 
as they are sequentially generated. Unfortunately, this procedure is poten- 
tially susceptible to the generation of a singular set of eigenvectors since 
it is very likely that no closed-loop matrix with the precise set of eigenvalues 
and eigenvectors as selected would exist for any choice of control law. This 
is exactly what occurred in example 2 of the section entitled "Numerical 
Examples." However, is it possible to detect the occurrence of such a situa- 
tion by selecting the eigenvectors sequentially as follows. 

Let N = ^vi : V2 : . • . : be the matrix of k-1 linearly 

independent eigenvectors already generated. Then the projector of the sub- 
space spanning these eigenvectors is given by 

R(k-l) = n[nTn] 


Selecting the kth eigenvector v^ so that 


R(k-l)v|^ ^ Vk 


(15) 


insures the linear independence of the kth eigenvector. In the numerical 
example 2, equation (15) cannot be satisfied for the third eigenvalue and eigen- 
vector specification. This condition implies a slight perturbation in the 
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eigenvalue and/or eigenvector specification must be made to insure the linear 
independence of the eigenvectors. 

An algorithm is now presented which incorporates the condition of equa- 
tion (15) without having to compute the projection matrices explicitly 

The sequential synthesis of the eigenvectors is accomplished in a special 
canonical form in which the testing of the linear independence between the 
eigenvectors is straightforward. Further, the procedure allows the designer 
to insure the linear independence between the synthesized eigenvectors to a 
degree determined by a prespecified tolerance parameter. For clarity of pre- 
sentation, the algorithm is detailed for real eigenvalues. The extension to 
complex conjugate pairs follows as a natural extension of the real eigenvalue 
problem. 

The following notations are used throughout the algorithm presentation: 

(1) e^ is an n-vector with rth entry equal to unity and all other 
entries zero. 

(2) v^T = [z^T . „^T] is the kth eigenvector, where Z]^ is the designer 
specified m-vector. 

(3) where q( 0) = and i / 0, will be 

defined in notation (6) . 

(4) is the rth element of and for convenience 

of notation 

(5) M(k) is an n X n matrix of the form 



In-r 

« • • 

. ma 

0 

* « • 



M(k) = 

0 

1 

0 


(16) 


0 


• • • 

Ir-1 



where the vectors 


and 

are computed so that 

the vector 

is transformed to 

a canonical form as 


' 


M(k)v^(k-1) = 

Oker,' r e 

(k)^ 



where is a 

subset 

of integers 

{1,2,.. 

.,n} containinc 

the indices not 


already used in the construction of the matrices . . ., 
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and A(l) is the complete set Note can be constructed, 

if and only if, 0]^ f 0. 

(6) Q(k-l) = M(k-l)M(k-2) . . . (k > 1) 

The algorithm now proceeds as follows: 

Step 1: For k = l,2,...,n, do steps 2 to 5. 

Step 2: For X = X]^, compute det -m “ f]* 

(a) If det = 0, perturb Xj^ to (Xj^ + 6X|^) and repeat step 2. 

(b) If det ^ 0, go to step 3. 

Step 3: For X = Xj^, compute C|^ (eq. (10)) as 

Ck = [^kln-m - + M] (17) 

Step 4: For some r e ^A(k)^ 

(a) Compute the m- vector 

[gr(k)]’’ = f,(k-l) + h^tk-DCk 

where ' hj-(k~l)) is the rth row of the transformation 

g(k-l) (notation (6)) and h^ (k~l) is a (n - m)-row vector. 

(b) Compute 

^k = Zk (18) 

where is the arbitrarily specified m design vector. 

(i) If Ok f 0, compute Wk (eq. (10)) and n(k) (eq. (16)) 
and go to step 1. 

(ii) If Ok = 0, select another r e ^A(k)^ and return to 
step 4 (a) . 

(iii) If Ok = 0 for all r e ^A(k)^, go to step 5. 

Step 5: For some r e (A(k)} 

(a) If gr (*^^ ^ 0, perturb Zk to (Zk + 5zk) to make ^k 5^ 8 
(eq. (18)), compute Wk and M(k), and go to step 1. 
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(b) If gjf = 0, select another r e and repeat step 5(a). 

(c) If = Of all r ^ perturb Xj^ to (Xj^ + SXj^) and 

go to step 2 . 

Step 6: Compute the feedback gains by using equations (12). 

Step 7: Stop. 

The kth linearly independent eigenvector V|^ can be synthesized pro- 
vided 0]^ ^ 0. This can be seen clearly as follows. Without loss of gen- 
erality, let the first (k - 1) eigenvectors be generated with indices 
r = 1,2,..., (k - 1). Then these vectors are transformed into a canonical 
form under Q(k“l) as 


Q(k-l) 

: V2 



Diag [cfi,02,.-.,0k-l] 
0 


with 


^ 0 


(i = l,2,»..,k — 1) 


and the projector spanning the subspace of these transformed eigenvectors has 
the simple form 


R(k-l) 



Now choosing the kth eigenvector so that its transformed vector becomes 
Vk^k-1) = Q(k-l)y|^ causes 

CTk ^ 0 (r ^ {k,k+l, . . . ,n}) 

and insures the linear independence of v^ since the constraint (eq. (15)) is 
clearly satisfied. 

The following observations can be made regarding the algorithm outlined; 

(1) The algorithm can be directly extended to complex pairs in quasi- 
diagonal form as described. Two eigenvectors are synthesized in one iteration 
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and equation (A5) is used instead of equation (8). Further, the test condi- 
tion 0)5 ^ 0 in step 4(b(i)) is modified, to testing the nonsingularity of 
a 2x2 matrix This matrix is constructed by selecting two rows of the 

transformation matrix obtained in the previous iteration and develop- 
ing a condition similar to equation (18). The transformation matrix for 

the complex pair is then obtained as a product of two transformations corre- 
sponding to the real eigenvectors : vi^+il associated with Pj^. This pro- 

cedure will insure the linear independence between Vj^ and vj^+i* 

(2) Explicit evaluation of the eigenvalues of matrix F is not needed to 

detect coincident mode assignment. It is sufficient to determine the appropri- 
ate determinant in step 2, this determinant being available as a byproduct in 
the synthesis of the kth eigenvector when C)^ is evaluated. If a mode is 
coincident with the spectrum of F, then a perturbation of the mode is required 
to insure [^X) 5 ljj_ m - F] is nonsingular and hence that equation (8) has a solu- 
tion for any arbitrary The degree of perturbation needed depends on the 

numerical tolerance set on the evaluation of the determinant. Also since the 
eigenvalue shift in step 2 is a designer's choice, system stability is always 
insured. Alternatively, if the eigenvalue assigned coincides with the spectrum 
of F, special eigenvector structures can be derived by noting that equa- 
tion (8) has a solution for = 0, and also for [|g + X^sjzi = 0 if 

(n - m) < m. (See ref. 6.) 

(3) The iterative procedure in step 4(b(ii)) attempts to meet exact 
eigenvalue/eigenvector specifications. In step 5(b) an attempt is made to meet 
exact eigenvalue specifications with slightly relaxed eigenvector specifica- 
tions (z-vector). The test in step 5(c) indicates that the eigenvalue specifi- 
cation implies that the corresponding eigenvector will lie in the eigensubspace 
already synthesized; this condition demands a perturbation in the eigenvalue 
specification. 

(4) Since the matrices have only one nontrivial column, coordinate 

transformations Q(k-l) reduce to simple vector multiplications. Further, the 
inverse of V required in equation (10) to evaluate the feedback gains in 
equations (12) is easily evaluated by noting that q(J^)v has the general form 


g(n)v = Diag [ai,02, . . . ,0n]L 


The matrix L is an elementary permutation matrix dependent on the 
sequence of generating the indices r e in steps 4 and 5, the 0 )^ are 

the nonzero pivotal elements, and 


Q(n) 


V"^ = L“1 Diag 


1 1 


'1 ^2 


'n 
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Also notice that 


det 



n 


TT 

k=l 


since 



and 


det [l] 


are unity. 


Thus, the numbers provide a good measure of the linear independence 

between the eigenvectors provided the eigenvector entries are scaled to a stan- 
dard basis as Vi^V£ = 1 for example. Thus, the numerical ill-conditioning 
of the modal matrix V for inversion is controlled by setting a tolerance on 
(k = l,2,...,n) to pass the test ^ 0. The important bearing of this 
mc^al matrix numerical condition on system sensitivity properties is discussed 
in the section entitled "Modal Sensitivity Analysis." 


(5) A noteworthy feature of the algorithm is that the eigenvectors do not 

explicitly undergo any change in the sequence of transformations Q(k-l). thus, 
the mode coupling characteristics of are kept transparent during synthesis, 

a very desirable feature for an off-line synthesis problem. 

(6) The algorithm outlined can be conveniently programed as an iterative, 
interactive multivariable synthesis procedure. 


Example 2 given in the section entitled "Numerical Examples" is now used 
to highlight the features of the spectral synthesis algorithm. Applying the 
algorithm step by step to the problem of example 2 yields the following 
synthesis sequence: 

(1) First mode (k = 1) : 


Xi = -1 

ziT = (1 0) 

A(l) = {1,2,3} 
Choose r = 1. Then 

Cl = (1 1) 

gi(l) = (1 0)T 

01 = 1 
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Since ^ 0, the first eigenvector can be synthesized as = (1 0 l)"^; 

and the transformation matrix q( 1) becomes 


1 


q(1) 


0 


-1 


0 0 
1 0 
0 1 


(2) Second mode (k = 2) ; 


X2 = -2 
Z 2 T = (-1 1) 

A(2) = {2,3} 
Choose r = 2. Then 

C2 = (0.5 0,5) 

g2<2) = (0 1)T 

02 = 1 


Since O2 f 0, the second eigenvector can be synthesized as 
and the transformation q( 2) becomes 


q(2) 



1 0 
1 0 


-1 


1 


1 


(3) Third mode (k = 3) : 


= — 1 


(t 

i' 

V2 = (-1 1 0)T; 


23 



(-1 1 ) 


Z3^ = 


A(3) = {3} 


Thus r = 3. Then 


C3 = (1 1) 

, 3 ( 3 ) =[0 0 ]'' 

CJ3 = 0 

Since 03 = 0, the third eigenvector cannot be synthesized. Further, since 
g 3 ( 3 ) = I'o 0]*^, the eigenvalue specification cannot be met. This result is 

obvious since the "basis" vectors spanning the two-dimensional subspace corre- 
sponding to = -1 are 

”1 0" 

0 1 

1 1 


and {vi,V 2 } already span this subspace. Thus, a second eigenvector associated 
with X = -1 cannot be synthesized as revealed by the null vector g 3 ^^^. This 
result implies X 3 must be perturbed slightly. Let X 3 = (-1 - e) ; e > 0. 

Then 


C 3 = 


>1 + e 1 + e. 


93 


(3) = 


-e 


-e 


,1 + e 1 + e, 


and again 
this case, 
can be met 
Z3T = (-1 


03 = 0; thus, the third eigenvector still cannot be synthesized. In 
since g 3 ^ 3 ) is not the null vector, the eigenvalue specification 
but the eigenvector specification requires slight perturbation. Let 
1 + 6 ) ; 6^0. Then 
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-e6 

03 = — 5^ 0 

1 + e 


and V3 can now be synthesized as 


V3 = -1 1+6 


1 + ei 


The closed-loop modal matrix V is 


1 -1 -1 


V = 


0 


1 


1 1+6 

6 

0 

1 + e 


and A = Diag [-1 -2 



Further 


det 



6e 

1 + e 


since = 02 = !• Thus, in this case, setting a tolerance on the value of 
03 would directly control the numerical ill-conditioning of V and conse- 
quently influence the choice of the perturbations e and 6. 


It is also interesting to note that if the sequence of assignment of modes 
were changed to = -1, X2 = - 1 , and X3 = -2, then both eigenvectors corre- 
sponding to Xi = -1 could be synthesized as 


1 . 




-1 

1 

0 
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and X3 = -2 could still be assigned without perturbation since 

= (-0.5 -0.5)'^. However# the eigenvector specification could not be 

met since 03 = 0. Thus, a perturbation in Z3^ would allow completion of 
the synthesis. 


MODAL SENSITIVITY ANALYSIS 

One of the primary objectives of a feedback control design is to insure 
that system responses remain well behaved under variations in plant parameters. 
It is possible to derive quantitative robustness measures for eigenvalue and 
eigenvector perturbations due to plant variations using the following sensi- 
tivity analysis detailed in reference 7. 

The closed-loop system satisfies 


Avi = X^vi (i = l,2,...,n) (19) 


Then for first-order differential changes in equation (19) , 

dA Vi + A dvi = Xi dvi + dXi Vi (20) 

For real eigenvalues let ti be a dual vector so that 

A^ti = Xiti (21) 

By premultiplying equation (20) by ti'*’ and using equation (21) the real eigen- 
value perturbation can be expressed as 


ti'*' dA Vi 

dXi = (22) 

tiTvi 


The corresponding eigenvector change is given by 


n-1 

dvi = ^ dijvi 

j=l 


(23) 
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with 


r 


(i = j) 


<3ij = < 




tjT(dA)Vi 
(Xi - Xj)(tjTvj) 


(i j) 


Similar expressions can be easily developed for complex conjugate pairs. (See 
ref. 6.) From equations (22) and (23), bounds can be established for the 
eigenvalue/eigenvector changes as 


dXil S 


dA 


^ill II ti 


ti^Vi 


= -^illdAl 


where 

Vi II ||t 

ti^Vi 


(24) 


(25) 


where is defined as the "mode condition number" corresponding to the eigen- 

value X^ and clearly 71^ ^ 1. if the eigenvectors V£ are normalized so that 

ViTvi = ti^vi = 1 (26) 

Then = ||ti ||2 and 


IdXil S Tri||dA||2 (27,) 

where itill 2 = ti'^ti and ||dA ||2 = ^Maximum eigenvalue of [dA'^ . dA^. 
Similarly, the corresponding bound for eigenvector perturbation is given by 


n 


||dvill2 ^ ||dA||2 ^ 


j=l 


TT 


j 
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The condition numbers irj[ in equation (27) provide quantitative robust- 
ness measures for the nominal eigenvalues and eigenvectors of the closed-loop 
system A. It is also clear that for plant perturbations not known a priori, 
the system has best robustness properties if = 1 (i = l,2,...,n), which 
implies that the nominal modal matrix V is orthonormal with t^ = V£ 
(i=l,2,«>«,n)> 

From equation (8) it is apparent that sufficient eigenvector freedom is 
generally unavailable to construct an orthonormal modal matrix. Thus, from a 
practical design perspective it is better to seek a mode-decoupled design, with 
each state variable dominantly displaying one mode (as determined by a dominant 
eigenvector entry). This design, in the limit when total decoupling^ is possi- 
ble, results in an orthonormal V. Further, if the perturbations dA are 
known a prigri, then the appropriate zero sensitivity eigenvectors span the null 
space of dA, and if such an eigenvector can be synthesized by using state ^ 
feedback, the corresponding eigenvalue will be invariant to the specified dA. 

The mode condition numbers iti derived as the norm of the dual vectors 
^i (eq. (27)) also directly relate to the numerical ill-conditioning of the 
modal matrix V for inversion (linear dependence of eigenvectors) . This fact 
becomes apparent by considering the relation 


tl^ 

t2'^ 


. . . [vi ; V2 : . . . : Vn] = In 


(28) 


t T 
‘-n 


and noting that as the eigenvectors v^ tend toward being linearly dependent, 
the norms of the corresponding dual vectors t^ increase in order to satisfy 
the constraint equation (28) . Thus, large values of 11^ indicate that the 
system has poor robustness properties and also reinforce the fact that the 
corresponding modal matrix is ill-conditioned. It should be recalled that the 
spectral synthesis algorithm described earlier allows the designer to effi- 
ciently control this system robustness property by setting a numerical toler- 
ance on the 0]^ of equation (18) to pass the test 0 in steps 4 and 5 

of the algorithm. Since the determinant of V is obtained as the product of 
the a|^ (k = l,2,...,n), it is clear that the design process effectively con- 
trols the ill-conditioning of V for inversion. 
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AIRCRAFT LATERAL CONTROL SYSTEM SYNTHESIS 


In this section, the lateral performance objectives of an aircraft are 
formulated as an eigenvalue/eigenvector assignment problem and the utility of 
the spectral synthesis algorithm is demonstrated by designing a feedback con- 
trol system for a fighter aircraft. Major qualitative design objectives for 
lateral aircraft dynamics are 

(a) Fast roll rate response with minimum overshoot 

(b) Good Dutch roll damping 

(c) Low sideslip and peak lateral acceleration in response to 

roll stick command, for good turn coordination 

(d) Low roll response to sideslip gust inputs. 


Performance Specification 

The quantitative performance specifications can be briefly summarized as 
follows (ref. 8) : 

(1) The frequency and damping of the Dutch roll mode shall satisfy 


Cd > 0.19 
(Und >1*0 


Cd^nd ^ 0.35 


(2) If the spiral mode is unstable, the time to double shall be greater 
than 20 seconds. 

(3) The roll-subsidence time constant Tj^ shall be less than 1.0 second. 

(4) After a rudder-pedals-f ree step aileron command, the ratio of sideslip 
to the parameter 0 shall be less than that specified below. The aileron 
command shall be held fixed until the bank angle has changed at least 90°. 

(a) Category A (rapid maneuvering phase) 

3/0 < 6° for adverse sideslip 
3/0 < 2° for proverse sideslip 
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(b) Category C (take off /landing phase) 
3/0 < 10° for adverse sideslip 
3/0 < 3° for proverse sideslip 
The parameter 0 is given by 


fj^t 

90° 


for category A and 


(t = 1.3) 


4>t 

30° 


(t = 1.0) 


for category C, where <|)t is the bank angle achieved at t seconds after the 
step input. 

(5) For a step aileron command held until the bank angle has changed at 
least 90°/ the roll rate at the first minimum following the first peak shall 
be at least 60 percent of the peak value. 

(6) For category A turn coordination/ the lateral acceleration at the pilot 
shall be less than O.lSg for a 60°/sec roll to a 60° bank. 

It should be recognized that these aircraft response characteristics must 
usually be achieved under feedback gain magnitude constraints arising because 
of sensor noise and other considerations. For the numerical aircraft example 
to be discussed in the section entitled "Fighter Aircraft Control Analysis/" 
all the feedback gain magnitudes were constrained to be about unity. The air- 
craft performance in the numerical example is evaluated by using the linearized 
model at selected flight conditions. 


Aircraft Model and Design Considerations 

The nonlinear equations of motion of the aircraft are used to generate 
linear perturbation models at various flight conditions. Accordingly, the 
state space representation, referenced to the stability axes, takes the form 


X = Ax + Bu 


y = Cx + Du^ 


(29) 
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where 





Lg 

0 


^5a 

^«r" 

A = 


Nr 

Ng 

0 

B = 

< 

N6r 



Yr 

^3 

gAo 



^5r 


_ 1 

0 

0 

0 


0 

0 



1 

0 

0 

0 


0 

0 


0 

0 

0 

0 


0 

0 

C = 

Wp 

Mj. 

Mg 

M<t) 

D = 


M6r 


_ 0 

0 

0 

1 


0 

0 


with = (p lij 3 4>) ; uT = (6a ; 6j.) ; = (p \j) ay (j)) . 

For the aircraft control problem/ from sensor considerations, feedback 
from 3 variable is difficult to implement. Hence, the control laws are 
derived by using the output vector y of equation (29), where the lateral 
acceleration ay sensor substitutes for the sideslip 3 sensor. The state 
feedback law of equation (2) is modified as finding a control law of the form 


u = Ky + Up 


(30) 


such that the closed-loop system meets desired eigenvalue/eigenvector speci- 
fications, with K the output feedback matrix and Up the external pilot 
input. The closed-loop system after applying feedback law (30) takes the form 


X = A + BKC X + BP u, 


y = [C + DKCjx + DP u, 

where P = - Kd] ^ and K = PK is the equivalent output feedback 

matrix obtained by setting D = 0. K exists provided P exists. Further, 
the feedback law (30) is still equivalent to the state feedback case since C 
is rank n for all flight conditions. Thus, for every state feedback law 
u = Kx + Up (eq. (2)) derived by using the spectral synthesis algorithm, 
there exists an equivalent output feedback matrix K (eq. (30) given by 
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K = k[c + Dk]“^ 


(32) 


Direct matrix manipulations show that the inverse of the matrix [c + DI^ 
in equation (32) exists if P exists. Hence, the theory developed for the 
state variable feedback problem is directly applicable to the aircraft lateral 
dynamics model described by equation (29) . 

Appendix B gives the numerical values for the (A,B,C,D) quadruple of equa- 
tion (29) for a typical fighter aircraft at selected flight conditions. An 
examination of these models over the complete flight envelope reveals that the 
following system constraints contribute critically to the evolution of a satis- 
factory closed-loop system. 

The large Lg derived in most of the flight conditions indicates that 
significant roll motions are induced for sideslip gust inputs. A feedback 
augmentation to reduce this effect yields the following control equations; 




with 


and 



1 


Mg 




Mg = Mg + MggkAg + M5j.ki^ 


with the circumflex C) denoting closed-loop magnitudes. 

It should be noted that if the gain magnitudes k;^ and k^ are con- 

Sy 3y 

I 

strained, substantial reduction in the Lg magnitude are possible only if the 
magnitude of the control derivatives ^ large. Further, 

since the Mg magnitude is dependent on the ay sensor location, a possibility 
exists to optimally locate the acceleration sensor from a control viewpoint. 

For satisfactory turn coordination, 3 ay excitation at the pilot 

station must be minimized for aileron input. This minimization implies reduc- 
tion of cross coupling between the roll dynamics and the sideslip/yaw rate. 

This coupling arises not only because of mode cross coupling in the system 
matrix A, but also because of the input mixing through the control matrix B. 
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Whereas feedback is effective only in suppressing the intervar iable cross cou- 
pling in A, some form of precompensation is necessary to minimize the input mix- 
ing of B. Thus, Lj, Lg, Np, and Yp must be reduced by feedback. Again, 
the feedback magnitude constraints determine the limit of reduction. Thus, if 
the turn coordination specification cannot be met by feedback alone, a feed- 
forward compensation, in the form of an aileron/rudder interconnect, must be 
introduced to increase the ratio The resulting composite control law 

takes the form ' 


u = Ky + Hup 


(33) 


where Up = (6^ • is the external pilot input and H is a 2 x 2 non- 

singular feedforward matrix representing aileron/rudder interconnect. The 
closed-loop system after applying feedback law (33) takes the form 


where 


X = [a + BKc]x + BPH Up^ 

r 

y = [C + DKCjx + DPH Up 
K and P are defined in equation (31) . 


(34) 


Eigenvalue/Eigenvector Assignment Formulation 

The aircraft handling characteristics postulated earlier can be summarized 
as lateral dynamics being composed of two weakly coupled subsystems. Roll rate 
and bank angle constitute the first subsystem and display predominantly the 
roll subsidence and spiral modes. The second subsystem is characterized by a 
well damped Dutch roll mode defining the yaw rate and sideslip motions. This 
results in the mode (eigenvalue) and associated response variable (eigenvector 
shape) assignments of table I. 


TABLE I.- DESIRED MODAL SPECIFICATIONS 


Mode 

Dominant 

response variable 

Roll subsidence 

Spiral 

Dutch roll 

Roll rate, p 
Bank angle, (|) 

Yaw rate and sideslip, il^,8 


33 




TABLE II.- DESIRABLE CLOSED- LOOP MODAL STRUCTURE 


Eigenvector 

components 

Roll subsidence 
mode 

Dutch roll 
mode 

Spiral 

mode 

P 

1 

0 

Z13 

Z14 

'i' 

0 

1 

1 

^24 

6 

wii 

W12 

W13 

W14 


W12 

W22 

W23 

W24 


It should be emphasized that the structure in table II is only one of 
several configurations that can be selected. In practice, alternative struc- 
tures must be tried to effectively trade off different (often conflicting) per- 
formance specifications. 

In this report, the rationale for the selection of the Zj[j components 
in table II was to meet the mode decoupled structure postulated in table I. 

This structure, in addition to meeting the handling quality specification, 
would result in robust feedback designs. 

The roll subsidence mode design vector selection is straightforward. The 
Dutch roll mode eigenvector design parameters were selected to keep only one 
parameter (Z13) available for manipulation in the design iteration cycles. The 
parameter Z13 directly influences the closed-loop Lg magnitude and is 
selected to effect the desired reduction in Lg magnitude. The spiral mode 
eigenvector components z^4 and Z24 must be selected so that the bank angle 
response is dominated by the spiral mode. 

The g and (j> components (w^j elements) of the eigenvectors in table II 
are computed by using constraint equation (8) or (A5) , after specifying the 
appropriate eigenvalues. The eigenvalues must be selected to meet quantitative 
speed of response and at the same time yield desirable eigenvector forms. The 
mode coupling matrices Ci (eq. (10)) clearly define this constraint relation- 
ship, and aid in the selection of the appropriate eigenvalues. 

It is helpful at this stage to clarify these design concepts by using a 
numerical example. Consider the aircraft model corresponding to flight condi- 
tion 1 (appendix B) . In particular, consider the synthesis of the spiral mode 
eigenvector. The spiral mode eigenvalue can lie in the range -0.01 to -0.1. 

The design equations take the form (eq. (10)). For \ = -0.01, 
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W 14 ' 


-25.52 

-7.16' 


zi4 

W24 


-100 

-6.02 


*24 


( 35 ) 


Judicious choice of z ^4 and Z 24 can now be made by examining the modal 
coupling relations in equation (35) to meet the requirement that W 24 be very 
large compared with W 14 , Z]^ 4 , and Z24* acceptable solution is obtained 

by selecting Z 24 = - 3 . 5 z]^ 4 ; as a result, the normalized eigenvector (with 
V 4 ^V 4 = 1 ) is obtained as 


^ 4 *^ = (Z ]^4 Z24 W ]^4 W 24 ) = (—0.012 0.044 0.006 0.998) 


(36) 


which has the desired bank angle dominant structure. Suppose, on the other 
hand, that if the spiral mode is chosen as X = 0 . 1 , the design equations 
become 


Wi4 


~-4.52 

-11.59 


zi4 

_W24_ 


-10 

-0.6 


Z24_ 


and thus the eigenvector structure would be worse compared with that in equa- 
tion (36). Thus, a trade-off between speed of response and minimal influence 
of the slow spiral mode on the response of \|) and g would finally determine 
the spiral mode in the range -0.01 to -0.1. This analysis also highlights the 
fact that mode selection is not only dictated by speed of response but also by 
the eigenvector structure it generates through the mode coupling matrices 
(eq. (10)). In the design of the spiral mode eigenvector, it sould be empha- 
sized that suppression of this slow mode from the 3 response is vital in 
limiting steady-state 3 excursions under the roll-stick command. 

The synthesis procedure up to this stage can be summarized as consisting 
of the following steps: 

Step 1.- For each mode i compute the eigenvector coupling matrix 
(eq. (10)) for the range of assignable values. 

Step 2.- Select the mode from step 1 which yields the eigenvector struc- 
ture closest to the desired one by appropriately choosing the z vectors as 
illustrated in the numerical example. 

Step 3.- Complete the synthesis by specifying all the modes and design 
vectors and applying the spectral synthesis algorithm. 

It is very likely that a design in step 3 may violate gain magnitude con- 
straints. The design parameters must then be iteratively modified to meet 


35 



these constraints. The direction of parameter change is found by studying the 
open loop modal structure and identifying the design specifications that unduly 
violate the physics of the basic plant. For example, the 7^3 parameter 
directly determines the reduction of the Lg derivative and consequently con- 
trols the magnitude of the feedback gains k^ and k^ . 

3y 3y 

FIGHTER AIRCRAFT CONTROL ANALYSIS 

Twenty flight conditions describing the operational flight envelope of the 
aircraft are shown in figure 1. The design points are represented as a plot of 
Mach number against angle of attack. The altitude parameter is identified with 
different symbols. Linearized state variable models describing the lateral 
dynamics at each of these flight conditions were used for the synthesis anal- 
ysis. The flight conditions of figure 1 were classified into three categories 
as shown in table III. The group numbers are arranged in decreasing order of 
control effectiveness as measured by the magnitudes of the control derivatives 
L 5 and L(S , with group III flight conditions posing the worst conditions for 
synthesis from a mode decoupling viewpoint. 


TABLE III.- FLIGHT CONDITION CLASSIFICATION 


Flight group 

Flight conditions* 

I 

1, 2, 3, 4, 6 , 9, 10, 12, 13, 14, 19, 20 

II 

5, 7, 8 , 11, 15, 16 

III 

17, 18 


*See figure 1. 


Feedback Design 

The basic modal matrix structure of table II was employed for the flight 
conditions in groups I and II. The design parameter z^2 ranged from -0.2 
to -1 for group I and -2 to -4 for group II to meet the gain constraint require- 
ments on the lateral acceleration sensor. Larger magnitudes of 2^3 implied 
increased roll response to 3 gust inputs. As mentioned earlier, significant 
changes in the basic aircraft characteristics could not be achieved for low 
altitude and low Mach number flight conditions, namely, flight conditions 17 
and 18. Low control derivative magnitudes coupled with feedback gain limita- 
tions contributed to the problem. For these cases, the iterative design pro- 
cess was initiated with the Dutch roll mode eigenvector structure corresponding 
to that of the free aircraft. Final trimming of the design parameters and 
Dutch roll mode damping were effected to give a reasonable improvement in 
response with acceptable feedback gains. 
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Precompensation Design 

As noted earlier, if the turn coordination specification cannot be achieved 
by mode decoupling alone (through feedback) , additional feedforward compensation 
by the aileron/rudder interconnect is required. In order to meet the roll stick 
input specification, the following structure for the precompensation matrix H 
(eq. (33)) was found to be satisfactory 


H = 



0 

1 


(38) 




The interconnect gain 


kp, is determined to maximize the ratio of 


CONTROL AUGMENTED AIRCRAFT PERFORMANCE 

The modal control procedure was used to obtain the gain matrices K and 
H (eq. (33)) at all 20 design flight conditions summarized in figure 1. The 
next objective is to evolve a control system which operates over the flight 
envelope based on these fixed point designs. This procedure usually takes the 
form of developing gain schedules using appropriate air data parameters as 
scheduling parameters. However, in this report the question of gain scheduling 
will not be pursued; instead, the improvement in system performance, based on 
fixed point designs is shown at three selected design flight conditions 1, 

17, cind 20. These flight conditions cover conprehensively the complete flight 
envelope from the synthesis complexity viewpoint. Flight condition 1 represents 
the nominal cruise condition with no significant synthesis constraints. Flight 
condition 20 is a high-angle-of-attack condition usually known to result in 
poor lateral control and turn coordination (ref. 8). Finally, flight condi- 
tion 17 covers the landing approach condition and represents the most diffi- 
cult flight regime from a mode decoupling viewpoint. 

Tables IV to VI detail the eigenvalue/eigenvector modifications achieved 
at these flight conditions. In tables IV to VII, j = The results are 

based on the control law given in the respective tables. It is worthwhile 
noting the modification achieved in the eigenvector pair corresponding to the 
Dutch roll mode in flight conditions 1 and 20. The desired decoupling of the 
roll rate variable from the Dutch roll mode could be achieved only by assign- 
ing an overdamped complex pair of eigenvalues (-1.5 ± jO.75) as the Dutch roll 
mode. Further, the design iterations revealed that the Dutch roll mode damp- 
ing could not be reduced (to improve 'I' and 3 responses) without violating 
feedback gain limits on the ay sensor. 
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TABLE IV.- MODAL CHARACTERISTICS (FLIGHT CONDITION 1) 


[i - 


Eigenvector 

components 

Characteristic for free aircraft 
eigenvalue of - 

Characteristic for augmented 
aircraft eigenvalue of - 

-3.70 

-0.03 

-0.34 ± j2.66 

-6.00 

-0.01 

-1.5 ± jO.75 

P 

-0.964 

-0.032 

-0.403 0.829 

0.986 

0.013 

0. -0.117 


-0.041 

0.044 

-0.096 -0.131 

0. 

-0.045 

0.586 0.586 

B 

0.002 

0.002 

0.069 -0.034 

-0.008 

-0.0004 

0.137 0.527 

<t> 

0.261 

0.998 

0.314 0.112 

-0.164 

-0.999 

-0.041 0.034 


CONTROL LAW (eq. (33)) 

-0.204 -0.491 -0.799 0.017 

-0.152 0.318 -0.504 -0.021 


H 


1 

0.325 


0 

1 


TABLE V.- MODAL CHARACTERISTICS (FLIGHT CONDITION 17) 

[j = >R] 


Eigenvector 

conponents 

Characteristic for free 
eigenvalue of - 

aircraft 

Characteristic for augmented 
aircraft eigenvalue of - 

-1.97 

-0.063 

-0.168 

± jl.594 

-4.5 

-0.07 

-0.5 ± jO.9 

P 

-0.89 

-0.078 

-0.63 

0.556 

0.976 

-0.086 

-0.49^ 0.49 


-0.039 

0.139 

-0.03 

-0.075 

0. 

0.140 

-0.074 -0.098 

B 

-0.001 

0.017 

0.122 

-0.005 

-0.017 

-0.004 

0.139 -0.173 

<t> 

0.454 

0.987 

0.381 

0.357 

-0.217 

0.986 

0.643 0.198 


CONTROL LAW (eq. (33)) 


K 


-1.07 


-0.179 


0.595 -1.259 -0.294 

0.362 -0.171 -0.084 


H 


1 0 

-0.069 1 


38 





TABLE VI.- MODAL CHARACTERISTICS (FLIGHT CONDITION 20) 


[j 


Eigenvector 

components 

Characteristic for free aircraft 
eigenvalue of - 

Characteristic for augmented 
aircraft eigenvalue of - 

1 

o 

• 

-0.096 

-1.034 ± j2.898 

-6.0 

-0.01 

-1.5 ± jO.75 


mm 

0.102 

-0.046 0.946 

0.986 


0. -0.116 



-0.022 

-0.022 0.005 

0. 


0.578 0.578 


-0.022 

-0.007 

0.077 -0.028 

-0.044 

0.0004 

0.113 0.544 

4> 

-0.801 

-0.995 

0.296 -0.087 

-0.164 

0.999 

-0.073 -0.066 


CONTROL LAW (eq. (33)) 


K 


-0.357 

-0.083 


-0.305 -0.847 -0.001 

0.581 0.051 -0.033 


H = 


1 

- 0.02 


0 

1 


Table VII depicts the improvement in the aircraft robustness character- 
istics in terms of the mode condition numbers introduced in the section entitled 
"Modal Sensitivity Analysis." Significant improvement in system sensitivity 
characteristics has been achieved in flight conditions 1 and 20. However, the 
improvement in sensitivity performance in flight condition 17 is marginal and 
is due to the inability to appreciably decouple Dutch roll mode from p and 4> 
variables. In table VII the determinant of the modal matrix V with eigen- 
vectors normalized as v^'^^v^ = 1 is also included. It is evident from 
table VII that the determinant of V reflects the aggregate effect of all 
the mode condition numbers (tt^) . The determinant value decreases as the mode 
condition numbers increase. 

Finally, time response histories for roll stick input and sideslip gust 
inputs are compared in figures 2 and 3. From the response curves it is appar- 
ent that control augmentation has reduced the cross coupling between the roll 
axis (p,4>) and the yaw axis (i|^,B). It is also important to note that in fig- 
ures 2(a) and 2(c), feedback augmentation has reduced the peak roll rate capa- 
bility for pilot aileron input. This reduction can be offset by introducing 
appropriate gain factors in the forward loop (h matrix (eq. (33))). It should 
be noted that the roll response to sideslip gust input has been significantly 
reduced in figures 3(a) and 3(c). Finally, for flight condition 20 the quan- 
titative B/0 performance specification could not be met but was significantly 
improved compared with the free aircraft performance. In all other flight con- 
ditions the feedback system met all the lateral handling criteria postulated 
earlier. 
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TABLE VII.- AIRCRAFT MODAL SENSITIVITY CHARACTERISTICS 


Flight 

condition 

Free aircraft 

Augmented aircraft 

Eigenvalue 

Mode condition 
number, Hi 

Eigenvalue 

Mode condition 
number , Hi 

1 

-3.70 

10.39 

-6.0 

1.065 


-0 . 03 

2.27 

-0.01 

1.027 


-0.34 ± j2.66 

14.09 

-1.5 ± jO.75 

4.337 

17 

-1.97 

6.76 

-4.5 

3.51 


-0.06 

3.52 

-0.07 

3.17 


-0.17 ± jl.59 

9.68 

-0.5 ± jO.9 

6.817 

20 

-0.77 

6.62 

-6.0 

1.069 


-0.096 

6.36 

-0.01 

1. 021 


-1.03 ± j2.9 
1 

14.32 

-1.5 ± jO.75 

4.0 


Flight 

Free aircraft 

Augmented aircraft 

condition 

Det [V] 

Det [v] 

1 

-0.0146 

-0.225 

17 

-0.0202 

0.0412 

20 

-0.0086 

0.242 


CONCLUDING REMARKS 

A new multivariable synthesis procedure for aircraft handling qualities 
design has been developed. The synthesis procedure is computationally simple 
in that it only involves solution of linear systems of equations. This results 
in an efficient computer-aided interactive design tool for flight control sys- 
tem synthesis. It is shown that the quantitative lateral design objectives can 
be conveniently formulated as an eigenvalue/eigenvector assignment problem. 

By noting the freedom available in the selection of eigenvalues and eigenvectors 
under state variable feedback, a systematic design methodology is evolved to 
meet the performance specifications. The inherently mode-oriented synthesis 
procedure provides significant insight into the design process and allows 
optimizing the response of the aircraft under often conflicting requirements. 

The study also revealed that the closed- loop modal matrix (matrix of eigen- 
vectors) plays a pivotal role in characterizing system sensitivity to plant 
parameter variations. New robustness measures related to the numerical ill- 
conditioning of the modal matrix for inversion, termed mode condition numbers, 
are also developed. 
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The synthesis technique is illustrated by using the linearized lateral 
dynamics model of a typical fighter aircraft. The simulation studies revealed 
that the feedback-augmented aircraft exhibited improved Dutch roll damping, 
good turn coordination for pilot roll stick command, and reduced roll response 
to sideslip gusts over the flight envelope. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
April 13, 1978 
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APPENDIX A 


COMPLEX CONJUGATE PAIR IN REAL ARITHMETIC FORM 

In order to perform real arithmetic computations with complex conjugate 
pair eigenvalues, the following transformation is established for reference. 

Consider the eigenvalue/eigenvector relation 


A(q + js) = (a + jb) (q + js) (Al) 


where a + jb is a complex eigenvalue of A and q + js the associated eigen- 
vector (with j = \pl) . Equation (Al) can be solved for real and imaginary 
parts as 


Aq = aq - bs" 


As = bq + as^ 


(A2) 


with the same relation holding for the conjugate eigenvalue a - jb. Equa- 
tion (A2) can be written as 


/a -b\ 

A(q : s) = (q : s)l 1 (A3) 


where q and s are real vectors corresponding to the real and imaginary part 
of the complex eigenvector and the complex conjugate pair of eigenvalues a ± jb 
can be written in equivalent real form as 


P 



b 

a 


(A4) 


By using the real arithmetic representation for a pair of complex eigenvalues 
®i - 3^i' the eigenvector constraints of equation (8) take- the form 
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APPENDIX A 







— ~ 

n-m “ ^ • “t>iln-m 

• 


Wi 


G + 

• 

zi 

i^n-m • ^i^n-m ” ^ 


wi+1 


• 

bj^S • G + 

Zi+1 


(A5) 


where 


[vi : Vi+i] 


Zi 


Zi+1 

Wi 


Wi+1 


are the real eigenvector pair associated with a^ ± jbi 
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APPENDIX B 


AIRCRAFT STATE VARIABLE MODEL AT SELECTED PLIGHT CONDITIONS 
The lateral dynamics of the aircraft modeled in state space takes the form 


X = Ax + Bu 

> (Bl) 

y = Cx + Du 

where x"^ = (p 4) 6 <j)) ; = (63 : 6 ^) ; and y = ay. For brevity of doc- 

umentation the state variables p, ij), and (J) are omitted from the output set 
in equation (Bl) . The respective quadruples (A, B, C, 6 ) for selected 
flight conditions are listed in tables Bl to B3. 


TABLE Bl.- FLIGHT CONDITION 1 (NOMINAL CRUISE) 

Mach number, 0.67; altitude, 6096 m (20 000 ft); 
angle of attack, 3.45° 


A matrix: 


-3.79E+00 

-1.34E-01 

6.02E-02 

l.OOE+00 


4.06E-02 

-3.59E-01 

-9.97E-01 

6.03E-02 


-5.20E+01 

4.24E-<-00 

-2.72E-01 

0 . 


0 . 

0 . 


4.62E-02 

0 . 


B matrix: 


2.50E+01 

1.42E+00 

5.01E-03 

0 . 


9.83E+00 

-4.20E+00 

5.03E-02 

0 . 


C matrix: 


-1.25E-01 


-6.12E-02 


-3.41E+00 


-1.50E-03 


D matrix: 

1.03E+00 -2.66E-01 
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APPENDIX B 


TABLE B2.- FLIGHT CONDITION 17 (LANDING APPROACH) 

l^ch number, 0.19; altitude, 30 m (98.43 f t) ^ 

|_ angle of attack, 6.72° J 


A matrix: 

-1.99E+00 

-7.97E-02 

1.19E-01 

l.OOE+00 


B matrix: 

2.98E+00 

2.47E-02 

3.50E-03 

0 . 


C matrix: 

-6.35E-02 5.93E-03 -9.79E-01 -4.38E-03 


D matrix: 

9.01E-02 -2.89E-02 


1.79E+00 

-7.92E-01 

3.74E-02 

0 . 


9.21E-01 -1.65E+02 0. 

-1.77E-01 5.68E-01 0. 

-9.92E-01 -2.01E-01 1.51E-01 

1.18E-01 0. 0. 
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APPENDIX B 


TABLE B3.- FLIGHT CONDITION 20 (HIGH ANGLE OF ATTACK) 

Mach number, 0.6; altitude, 6096 m (20 000 ft); load 
factor, g units, 4; angle of attack, 15.44° 


A matrix: 

-2.25E+00 

-7.03E-02 

2.66E-01 

l.OOE+00 


B matr ix : 

1.61E+01 

7.25E-01 

1.90E-04 

0 . 


C matrix: 

-1.04E-01 3.13E-01 -3.95E+00 -2.07E-02 


D matrix: 

5.76E-01 -1.93E-01 


8.21E+00 

-3.43E+00 

4.76E-02 

0 . 


7.58E-01 -3.63E+01 0. 

-4.06E-01 -4.50E-02 0. 

-9.63E-01 -2.83E-01 4.98E-02 

2.76E-01 0. 0. 
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Time, sec 

(a) Flight condition 1 (nominal cruise) . 

Figure 2.- Comparison of free aircraft and augmented aircraft response to roll 
stick step input. 6^(0) = 1°; A indicates free aircraft response and B cor- 
responds to augmented aircraft response. 
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Rudder deflection. Aileron deflection. Lateral acceleration, 
deg deg g units 




Figure 2.- Continued. 
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Rudder deflection, Aileron deflection, Lateral acceleration) 





(c) Concluded. 
Figure 2.- Concluded. 
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Time, sec 


(a) Concluded. 


Figure 3.- Continued. 







Time, sec 


(b) Concluded. 


Figure 3.- Continued. 
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